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Abstract 

This paper derives a new semi closed-form approximation formula for pricing an up- 
and-out barrier option under a certain type of stochastic volatility model including SABR 
model by applying a rigorous asymptotic expansion method developed by Kato, Taka- 
hashi and Yamada [TJ. We also demonstrate the validity of our approximation method 
through numerical examples. 
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1 Introduction 

Numerical computation schemes for pricing barrier options have been a topic of great 
interest in mathematical finance and stochastic analysis. One of the tractable approaches 
for evaluation of barrier options is to derive an analytical approximation. However, from 
the mathematical viewpoint, deriving an approximation formula by applying stochastic 
analysis is not an easy task since the Malliavin calculus approach as in Takahashi and 
Yamada [3j cannot be directly applied. Recently, Kato, Takahashi and Yamada pQ 
has provided a new asymptotic expansion method for the Cauchy-Dirichlet problem by 
developing a rigorous perturbation scheme in a partial differential equation (PDE), and 
as an example, derived an approximation formula for a down-and-out call option price 
under a stochastic volatility model. In this paper, we give a new asymptotic expansion 
formula for an up-and-out call option price under a stochastic volatility model which is 
widely used in trading practice. Moreover, we show the validity of our formula through 
numerical experiments. 
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2 Asymptotic expansion formula for up-and-out 
barrier option prices 

Consider the following stochastic differential equation (SDE) in a stochastic volatility 
model: 

dSf = (c-q)S!dt + a!SidBl 
S £ o = S, 

da\ = e\{6-o £ t )dt 

+eua £ t (pdB} + y/l - p 2 dB 2 ), 
erg = cr, 

where S,a,c,q > 0, e G [0,1), X,6,v > 0, p G [-1,1] and B = (B 1 ^ 2 ) is a two 
dimensional standard Brownian motion. This model is motivated by pricing currency 
options. In this case, c and q represent a domestic interest rate and a foreign interest rate, 
respectively. The process S £ denotes a price of the underlying currency. Our purpose 
is to evaluate an up-and-out barrier option with time-to- maturity T — t and the upper 
barrier price H(> S), and its initial value is represented under a risk- neutral probability 
measure as follows: 

^Barrier ~~ *> 



E 



" C(T /( 5 'T-i) 1 {T(o 1 H)(5 e )>T-t} 



where / stands for a call option payoff function f(s) = max{s — K, 0} for some K > 0. 
Here, the stopping time t^ ^(S £ ) is defined as 

r (0 ,H)(S £ ) = mf{t G [0,T};S! $ (0,H)} (inf := oo). 

Remark that CBarrier(-^ — *i nas no closed-form solution and therefore we have to 
rely on some numerical method such as the Monte-Carlo simulation in order to calculate 
Cg^r icr (T — t,S). However, when e = 0, C^® icT (T — t, S) corresponds to the up-and-out 
barrier option price in the Black-Scholes model which is known to be solved explicitly. 
Then, for e > 0, we are able to derive a semi closed-form expansion around Cg^ ier (T — 
t, S) when £ | 0. This is our main result and hereafter we show our approximation 
method for Cf^ r e icr (T - t, S). 

Clearly, applying Ito's formula, we can derive the SDE of logarithmic process of Sf 

as 

dXf = { c -q- l -{a £ t f)dt + a £ t dBl 
X £ = x := log S. 
Then we can rewrite Cg^ ier (T — t, S) as 

^Barrier (T — *' &X ) 

where f(x) = max{e x — K, 0} and D = (— oo, log Note that 
r D (X £ ) = m£{t G [0,T] ; Xf £ D} = r^ H) (S £ ). 
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Let u 6 (t, x) = Cgarrier^ ~~ eX ) ^ or * e IP> -Tl anc ^ ^ € R. Then u e {t,x) satisfies the 



following PDE: 



where 



d \ 

— +^ £ -cju £ (t,x)=0, (t, x) G (0, T] x D, 

u £ (T,x) = f(x), x & D, 

{ u £ (t,logtf) = 0, t£[0,T], 



( 1 i\ d 1 2 d 2 

2^ ... ,3 2 1 2 2 <9 2 



<9xcV 9<7 2 <9<7 2 

As mentioned above, when e = 0, we can obtain the explicit value of u°(t,x). In this 
case, u°(t,x) = C#f TTicT (T — t, e x ,a, H) represents the price of the up-and-out barrier call 
option under the Black-Scholes model. We have 



where 



with 



y 
yi 

See Hull [2] for the details. 



n BS _ riBS n 

'-'Barrier — ° Vanilla ~~ U ) 



C^ma = e x e-i T N(d 1 )-Ke- cT N(d 2 ), 
C = e x e~ qT N{x 1 )-Ke~ cT N{x 2 ) 

-e x e-^(^) 2 [N(-y)-N(- yi )\ 

x [N{-y + aVf) - N(- yi + aVf)} 

x - log H + (c - g)T + l/2<7 2 r 

xi — aVT, 

(c-q) 1 
A " ~^^ + 2' 

2 log if - x - log if + (c - g)T + 1 /2cr 2 T 



xi 

X2 



log H - x + (c - g)T + l/2o- 2 T 



We can represent u°(t,x) = Pf ) f(x) by using a semi-group (Pf)t defined as 

2(logH-g)(logH-B] 



P s D g(x) = I e— (1-e- 



t 

losH 
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1 <.y-x~(c~q~^cr 2 )s) 2 



--e ' g(y)dy 



(2.1) 

for a continuous function g with polynomial growth rate which satisfies g(x) = on dD. 

The main result of Kato, Takahashi and Yamada pQ suggests the following approxi- 
mation formula. 

[Asymptotic expansion formula] 



u £ (t,x) = C ] BS 



Barrier 

T-t 



where 



+ee-« T -* [ ~ P?^P?_ t _J{x)ds + 0{e% 
Jo 



Using (|2.ip . the term / P s D ^f 1 °P^?„ i „ s /(x)(is is expressed as follows: 
Jo 

j T ~ l ' P?^P°_ t _J{x)ds 
Jo 



T-t /-log H 2(lo e H-x) {log H-y) 

e" cs (l - e ^ ; 



J-oo 

(„ — !r.-l, n-i^ 2 Vs1 2 

3D 



I (y-x-(c-q-^CT 2 )s) 2 



(2.2) 



We are able to compute the integrand of the right hand side of the above formula 
([£2D as 

d d 2 
Here, —C^ VTlcy .(T,e x ) and - - C^„ icr (T,e x ) are concretely expressed as follows: 

^BarrierV- 1 ) ° / 



e- 9T e x n(di)v / r 

-e-« T e a: n(xi)v / r - (H - K) e - cT n(x 2 ) 



-Xl 
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x I (log H - x) ~ A{c ~ q) {N{-y) - N(- yi )} 



2A 



+[n(y)— - n(yi)— ] 
a a 



\e x 



x < (log H — x)- 



2X-2 



-4(c 



O"' 



.[N(-y')-N(-y[)] 



+ [n(y')- - n(y[) — 
a o 



d 2 
dxda 

-e~ qT e x n(xi){-X2) — 
a 

-{H-K)e-^^={ Xl x 2 -l} 

+ 4(c-g) + H — x) + 1} 

a 6 

+e X e- qT (*) 2X ny) y --n( yi ) y -±] 
\e x J a a 

-Ms)") 



ct 3 V VcVr/ \aVT 

2A 



x \ n{y)-^-j={yy' - 1) - rc(yi)-^=(yiyi - 1) 
-Ke- cT [iV(y')-^(y;)] 

|) 2A " 2 fc^{(2A-2)(log^-x) + l} 

+ Ke-^(^\lo g H-x^ C - q) 



\e x J a 
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x = -n{y[) 



+Ke~ cT (2X - 2 



jj\ 2A-2 



x \n(y')- - n(y[) — ) 
\ a a J 



-Ke~ cT [ - ) 



H\ 2A ~ 2 



x ( n{y')^-={y'y - 1) - \=.(y' x y x - 1) 



where 



1 



. 21ogif-x-logK + (c-g)T- icr 2 r 

y = 



, _ log # - x + (c - q)T - \a 2 T 

3 Numerical Examples 

In this section we show numerical examples for pricing European up-and-out barrier 
call options under SABR volatility model (A = 0) as an illustrative purpose. By the 
asymptotic expansion formula in the previous section, we see 



^Barrier C^"' ^) — ^Barrier 

+ee- cT f T P s D J%Pj>_J(S)ds. 
J o 

Let us define AE first and AE zeroth as 

AE first = Ci^ ei (T,3) 

+ee- cT [ T P s D J??P°_ s f(S)ds, 
Jo 

AE zeroth = Cf| rier (T, S). 

Below we list the numerical examples, [Case 1] - [Case 6], where the numbers in the 
parentheses show the error rates (%) relative to the benchmark prices of C-qJ^ CT, S) 
which are computed by Monte-Carlo simulations with 100, 000 time steps and 1, 000, 000 
trials (denoted by MC). We check the accuracy of our approximation formula by chang- 
ing the model parameters. 

Apparently, our approximation formula AE first improves the accuracy for (T, S), 

and it is observed that the approximation term ee~ cT l P® JC® Pt-sJ (S)ds accurately 

Jo 

compensates for the difference between C Bar '^ ier (T, S") and C^ rrier (T, S), which confirms 
the validity of our method. 

For all cases, we set S = 100, a = 0.2, c = 0.0, q = 0.0, p = -0.5, eA = 0.0, 6 = 0.0 
and T = 1.0. In Case 1, 2 and 3, given ev = 0.1, the upper bound price is set as 
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H = 120, 130, 140, respectively, while in Case 4, 5 and 6, given eu = 0.2, H is set as 
120, 130, 140, respectively. Particularly, for the case of eu = 0.2 (that is, higher volatility 
of volatility case, Case 4, 5 and 6), we remark that the errors of the approximation 
become slightly larger. However, as observed in comparison between AE first and AE 
zeroth, we are convinced that the higher order expansion improves the approximation 
further, which will be investigated in our next research. 

[Casel] 

S = 100, a = 0.2, c = 0.0, q = 0.0, eu = 0.1, 
p = -0.5, eX = 0.0, 9 = 0.0, H = 120, T = 1.0. 



Table 1: up-and-out barrier option prices and the relative errors (Case 1) 



Strike: K 


MC 


AE first 


AE zeroth 


100 


1.204 


1.188 (-1.35%) 


1.105 (-8.25%) 


102 


0.882 


0.869 (-1.44%) 


0.804 (-8.78%) 


105 


0.512 


0.504 (-1.62%) 


0.463 (-9.59%) 



[Case2] 



S = 100, cr = 0.2, c = 0.0, q = 0.0, eu = 0.1, 
p = -0.5, eX = 0.0, 9 = 0.0, H = 130, T = 1.0. 



Table 2: up-and-out barrier option prices and the relative errors (Case 2) 



Strike: K 


MC 


AE first 


AE zeroth 


100 


3.216 


3.200 (-0.49%) 


2.966 (-7.78%) 


102 


2.621 


2.607 (-0.55%) 


2.406 (-8.22%) 


105 


1.869 


1.857 (-0.69%) 


1.702 (-8.93%) 



[Case3] 



S = 100, a = 0.2, c = 0.0, q = 0.0, eu = 0.1, 
p = -0.5, eX = 0.0, 9 = 0.0, H = 140, T = 1.0. 
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Table 3: up-and-out barrier option prices and the relative errors (Case 3) 



Strike: K 


MC 


AE first 


AE zeroth 


100 


5.184 


5.186 (0.05%) 


4.847 (-6.49%) 


102 


4.420 


4.423 (0.06%) 


4.121 (-6.77%) 


105 


3.420 


3.422 (0.06%) 


3.174 (-7.19%) 



[Case4] 

S = 100, a = 0.2, c = 0.0, q = 0.0, ev = 0.2, 
p = -0.5, eX = 0.0, 9 = 0.0, H = 120, T = 1.0. 



Table 4: up-and-out barrier option prices and the relative errors (Case 4) 



Strike: K 


MC 


AE first 


AE zeroth 


100 
102 
105 


1.317 
0.971 
0.569 


1.271 (-3.51%) 
0.934 (-3.83%) 
0.545 (-4.30%) 


1.105 (-16.12%) 
0.804 (-17.15%) 
0.463 (-18.65%) 



[Case5] 

S = 100, a = 0.2, c = 0.0, q = 0.0, ev = 0.2, 
p = -0.5, eX = 0.0, 9 = 0.0, H = 130, T = 1.0. 



Table 5: up-and-out barrier option prices and the relative errors (Case 5) 



Strike: K 


MC 


AE first 


AE zeroth 


100 


3.475 


3.435 (-1.15%) 


2.966 (-14.66%) 


102 


2.844 


2.808 (-1.27%) 


2.406 (-15.42%) 


105 


2.041 


2.011 (-1.48%) 


1.702 (-16.58%) 



s 



[Case6] 



S = 100, cr = 0.2, c = 0.0, q = 0.0, ev = 0.2, 
p = -0.5, eX = 0.0, 9 = 0.0, H = 140, T = 1.0. 



Table 6: up-and-out barrier option prices and the relative errors (Case 6) 



Strike: K 


MC 


AE first 


AE zeroth 


100 


5.483 


5.526 (0.78%) 


4.847 (-11.59%) 


102 


4.683 


4.725 (0.85%) 


4.121 (-12.03%) 


105 


3.635 


3.670 (0.97%) 


3.174 (-12.68%) 
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